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ABSTRACT 

We present an analytical determination of the star formation rate (SFR) in molecular clouds, based 
on a time-dependent extension of our analytical theory of the stellar initial mass function (IMF). The 
theory yields SFR's in good agreement with observations, suggesting that turbulence is the dominant, 
initial process responsible for star formation. In contrast to previous SFR theories, the present one 
does not invoke an ad-hoc density threshold for star formation; instead, the SFR continuously increases 
with gas density, naturally yielding two different characteristic regimes, thus two different slopes in the 
SFR vs gas density relationship, in agreement with observational determinations. Besides the complete 
SFR derivation, we also provide a simplified expression, which reproduces reasonably well the complete 
calculations and can easily be used for quick determinations of SFR's in cloud environments. A key 
property at the heart of both our complete and simplified theory is that the SFR involves a density- 
dependent dynamical time, characteristic of each collapsing (prestellar) overdense region in the cloud, 
instead of one single mean or critical freefall timescale. Unfortunately, the SFR also depends on 
some ill-determined parameters, such as the core-to-star mass conversion efficiency and the crossing 
timescale. Although we provide estimates for these parameters, their uncertainty hampers a precise 
quantitative determination of the SFR, within less than a factor of a few. 
Subject headings: stars: formation — ISM: clouds — physical processes: turbulence 



1. INTRODUCTION 

The determination of the star formation rate (SFR) in 
molecular clouds and in galaxies is one of the main chal- 
lenges of star formation theory. In the modern paradigm 
of star formation, stars form out of prestellar cores which 
result from the gravo-turbulent fragmentation of molec- 
ular clouds (e.g. MacLow & Klessen 2004). Within the 
past few years, two analytical approaches have emerged, 
aiming at characterizing the SFR issued from the proba- 
bility density function (PDF) of density fluctuations in- 
duced by turbulence in the cloud (Krumholz & McKee 
2005-KM, Padoan & Nordlund 2011-PN). Both theories 
rely on (i) a density threshold, whose nature differs in 
each case, for star formation, (ii) one characteristic dy- 
namical timescale, defined either at the cloud's mean 
density or at the threshold density. In this Letter, we 
derive a SFR, based on our IMF analytical theory (Hen- 
nebelle & Chabrier 2008-HC08, 2009-HC09, Chabrier & 
Hennebelle 2011), and show that (i) this theory yields 
SFR values in good agreement with observations, (ii) 
there is no a priori density threshold for star formation; 
instead, the SFR continuously increases with gas density, 
with indeed two different regimes. We also show that the 
exact value of the SFR depends on the combination of 
some ill-determined parameters, notably the core-to-star 
efficiency and the crossing timescale, whose uncertain- 
ties, and dependence upon cloud conditions, hamper an 
exact determination of the SFR. 



2. STAR FORMATION RATE: THEORIES 

We first summarize the previous SFR theories by KM 
and PN. We then briefly present the SFR derived from 
a time-dependent extension of our theory of star forma- 
tion, which will be presented in details in a forthcoming 
paper. Finally we present a simplified version of this the- 
ory which, alternatively, can be seen as an improved KM 
or PN theory. 

Following Krumholz & McKee (2005), we define the 
dimensionless star formation rate per free-fall time, 
SFR ff , as the fraction of cloud mass converted into stars 

per cloud mean free-fall time, t^, i.e.: SFRff — ^f- rjL, 

where M* denotes the total star formation rate arising 
from a cloud of mass M c , size L c and mean density po- 

2.1. The Krumholz and McKee theory 

According to various simulations of hydrodynamic or 
MHD supersonic turbulence, the density PDF is well rep- 
resented in both cases by a lognormal form, 
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where M. is the Mach number and b ~ 0.5 (Federrath et 
al. 2010). 

The essence of the KM analysis is to assume that there 
is a critical density, p cr u, above which star formation is 
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occuring. Then, the SFR (eqn.(20) of KM) is simply 
obtained by estimating the fraction of gas with density 
larger than p crlU 



SF R ff = 
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with p = p/pa- KM further assume that r//. cr ~ rjf. 

In this expression, e is the (supposedly mass- 
independent) efficiency with which the mass within the 
collapsing prestellar cores is converted into stars. Cal- 
culations (e.g. Matzner & McKee 2000, Ciardi & Hen- 
nebelle 2010) as well as observations (e.g. Andre et al. 
2010) suggest that e ~ 0.3 — 0.5. The parameter <j> t corre- 
sponds to the time needed for a self-gravitating fluctua- 
tion to be replenished. KM estimate it to be of the order 
of a few, in agreement with the analysis we propose in 
the appendix. 

In KM, p cr a is determined from the condition that 
the corresponding Jeans length must be equal to the 
sonic length. Their underlying assumption is that tur- 
bulent support will be too efficient to enable star forma- 
tion at scales larger than the sonic length. This yields 
pcrit.KM = (0xA,/o/A s ) 2 , where 4> x is a coefficient of order 
unity, Ajo is the Jeans length at the mean cloud density 
and A s is the sonic length. 

2.2. The Padoan and Nordlund theory 

The expression obtained by PN is similar to the KM 
one, stated by eqn. fl2J, except that they do not assume 

r/Tj f = yjpcnt , as indeed 



T ff = T //> cr ' ms tead Tff jCr j 

comes out from the integral in eqn. @. They consider 
that both e and <j>t are equal to 1, except in the magne- 
tized case where they argue that e ~ 0.5 (which appears 
to be the main reason for the reduced SFR in the magne- 
tized case). With eqns. ((TJ and ^2j) , this yields (eqn. (30) 
of PN) 
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The main difference with the KM model, however, re- 
sides in the choice of pcrit- In PN, this latter is ob- 
tained by requiring that the corresponding Jeans length 
be equal to the typical thickness of the shocked layer, 
inferred by combining isothermal shock jump conditions 
and a turbulent velocity scaling v cx I 5 . This yields 
Peru, pn — 0.067 8~ 2 a V i r M 2 , where 9 m 0.35 is the ra- 
tio of the cloud size over the turbulent integral scale 
and a V i r is the virial parameter, a V i r = 2E^- m / E grav = 
5Vq 2 / (nGpoL 2 ), where Vq is the rms velocity within the 
cloud, representative of the level of turbulent vs gravita- 
tional energy in the cloud (eqns. (8-9) of PN). 

2.3. The Hennebelle and Chabrier theory 

In the HC theory of star formation (see HC08, HC09) 
prestellar cores are the outcome of initial density fluc- 
tuations that isolate themselves from the surrounding 
medium under the action of gravity. These fluctuations 
are determined by identifying in the cloud's random field 
of density fluctuations the structures of mass M which 
at scale R are gravitationally unstable, according to the 
virial theorem. This condition defines a scale- dependent 



(log)-density threshold, S R = \n(p c (R)/ p ), or equiva- 
lently a scale-dependent Jeans mass, M R 
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where C s is the sound speed, G the gravitational con- 
stant, aj a constant of order unity while Vq and r/ ~ 0.4 
determine the rms velocity: 
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A fluctuation of scale R will be replenished within a typi- 
cal crossing time tr, and will thus be replenished a num- 
ber of time equal to Tf f /(^ t rR,ff), where tr,// = m/fa 
is the freefall time at scale R (see appendix) , i.e. at den- 
sity pn ~ Mr/R 3 . Including this condition into the HC 
formalism yields, after some algebra, for the number- 
density mass spectrum of gravitationally bound struc- 
tures, M{M) = d(N/V)/dM: 
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Apart for the time ratio this equation is similar to 
equation (33) of HC08 and^equation (27) of HC09. Ac- 
cording to this definition, SFRff is thus given by the 
integral of the mass spectrum specified by equation ([6]): 
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According to eqn. (j4j, M cu t correponds to the mass as- 
sociated with the largest size fluctuations that can turn 
unstable in the cloud, y cu t = 2R/L C . We verified that, 
as long as y cu t is not too small, the results depend only 
weakly on its value (Hennebelle & Chabrier, in prep.). 
In the following, we will pick y cut ~ 0.1 as our fiducial 
value. 

A few remarks are worth discussing at this stage. First, 
in the present theory there is no explicitly introduced 
critical scale or density for star formation, as we sum 
up over all gravitationally unstable cores, irrespectively 
of their scale or density. This is achieved through the 
multiscale analysis expressed by eqns. (j6])-<l7j) . Indeed, 
turbulence is by essence a multi-scale phenomenon and 
introducing a critical scale does not appear clearly jus- 
tified. Indeed, a piece of fluid, even if dominated by 
turbulence, can still collapse if it is self-gravitating. 

Another essential difference with the KM and PN the- 
ories is that these latter rely on a unique characteris- 
tic collapsing time (the mean cloud freefall time in KM 
and the critical density freefall time in PN) . This can be 
seen from eqns. ([2]) and ([3]), where the term p 1 ^ 2 , taken 
at Pcrit i lies outside the integral. In constrast, in our 
theory (see eqn. ([7])), the freefall density dependence of 
each collapsing structure is properly accounted for, as the 
freefall time consistently varies with mass M and scale 

R, tr ff oc p R . 
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2.4. A simplified multi-freefall theory 

Even though we stress that the SFR cannot be properly 
determined by a simple integral of the density PDF, be- 
cause such an integral, unlike eqn. ([7]), does not take into 
account the spatial distribution of the gas, we suggest 
the following simplified but more consistent expression, 
which retains the collapsing time density-dependence, in- 
stead of eqns. ((2]) and (j3|): 
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This SFR is larger than the ones given by eqns. © and 
([3]), as shown in the next section. 
We consider different choices for p cr it- When using 

Pcrit, KM 01" Pcrit 

t PN, we refer to the corresponding SFR 
as " multi-freefall KM or PN" , respectively. We also con- 
sider another value for p cr i t , obtained by simply requiring 
that the Jeans length at this density is equal to y cu tL c - 
This corresponds to the assumption that only fluctua- 
tions smaller than a given cloud size fraction can col- 
lapse. We simply refer to this model as "multi-freefall." 
It is easy to check that eqn. ((8]) depends only weakly on 
y cuU except when y cut -» 0. 

3. RESULTS 

We now compare the SFR / / predicted by the various 
theories and confront the results to recent observations. 

3.1. Comparison between the various theories 

As in KM and PN, we define the cloud properties by 
a V i r and M.. Figure CD displays SFR° ff , which corre- 
sponds to SFRff for e = 1 and (fit = 1, obtained with 
different formalisms, for various values of the virial pa- 
rameter, and for three typical Mach numbers, namely 
M =16, 9 and 4. 

Both the HC and the multi-freefall models are larger by 
a factor ^2-3 than PN and by at least an order of mag- 
nitude than KM. This stems from the fact that, when 
taking into account the density-dependence of the struc- 
ture collapsing times, (i) dense regions collapse fast, and 
(ii) fluctuations denser than p cr u have a smaller free-fall 
time than the ones at p cr u , globally increasing the value 
of SFRjj:. When such a density-dependence is prop- 
erly accounted for, all SFR determinations are in bet- 
ter agreement. Nevertheless, some differences persist be- 
tween the various multi-freefall models. This stems from 
the choice of p cr it- Indeed, the choice of p cr it in the KM 
and PN theories (see §2.1 and 2.2) yields y cut ~ M~ 2 
and thus corresponds to very small values of y cu t , imply- 
ing that only small Jeans masses (or conversely only very 
dense structures) are taken into account in these models. 

Two interesting trends can be inferred from Fig. 1. 
First, increasing the virial parameter leads to a decrease 
of the SFR, with a severe reduction above some typi- 
cal value of a v i r , which decreases with decreasing Mach 
number. This naturally arises from the fact that, 
increases, the increasing contribution of kinetic energy 
over potential energy prevents gravitational collapse and 
thus inhibits star formation, a point already noticed by 



KM and PN. Second, the SFR increases, although mod- 
estly, with the Mach number. This is because increasing 
the Mach number extends the core mass function (CMF) 
into the low-mass domain (see PN and HC08) and, as 
small-scale structures have shorter free-fall times, this 
increases the number of collapsing small cores, thus the 
SFR. This positive dependence of the SFR upon M. is 
in agreement with the results of PN but contrasts with 
the ones of KM (see their equation 30), as seen in the 
figure. Such a decreasing dependence of the SFR with 
increasing Mach in the KM theory clearly stems from the 

fact that the p crit term (see eqn. @) is lacking in their 
eqn. (20). 

3.2. Comparison with observations. 

Star forming giant molecular clouds in the Milky Way 
have masses 10 3 <M c /M o <3 x 10 6 , with M ps 4-30 
and a V i r ps 0.3-3. The observed SFR per cloud free-fall 
time lies in the range 0.03 < SFRf f < 0.3, with a mean 
value (SFR ff ) ps 0.16 (Murray 2011, although see Feld- 
mann & Gnedin 2011 for caution). Evans et al. (2009) 
and Heiderman et al. (2010) find SFR's in the range 
ps 0.02-0.12 for nearby molecular clouds and ps 0.03-0.5 
for massive star-forming dense clumps, yielding a mean 
value ps 0.1, about an order of magnitude larger than the 
values predicted by KM. Krumholz & Tan (2007, Fig. 5) 
report lower values at low density. At high density 10 
cm -3 ), however, two of their three data (ONC and CS(5- 
4)), are compatible with the aforementioned mean val- 
uefl According to our calculations (see Fig. 1), for such 
cloud/clump characteristics, SFRff is predicted to lie 
within the range ps 0.3-3. As discussed in § 12. H and in the 
appendix, the effective SFR is SFR ff = (e/<fi) x SFR° ff , 
where e/cfit ~ 0.1-0.2. Therefore, according to the present 
calculations, theories based on a multi-freefall formal- 
ism yield SFR's per free-fall time in typical molecular 
clumps in the range SFRff ps 0.03-0.6, going from low- 
dense clouds to the densest clumps, consistent with the 
observed values. 

Figure 2 displays the SFR per unit area, = 
SFRff x Ytg/r^f, with S g = L c p = M c /ttL 2 c , as a func- 
tion of cloud surface densities, E g , for four typical cloud 
sizes, L c = 1, 4 10 and 40 pc. The clouds are assumed to 
follow Larson's (1981) relations and thus have velocities 
given by eqn. ([5]) and densities no x (-^c/lpc) -0 ' 7 , where 
n = 10 2 to 10 4 cm -3 , yielding cloud masses M c = 200- 
10 6 Mq. From these values, M. and a V i r can be consis- 
tently determined. In Fig. 2, we have taken e/cfit = 0.10 

Lada et al. (2010) and Heiderman et al. (2010) data 
for clouds and massive clumps are shown for comparison. 
Clearly, almost all theories exhibit a direct correlation 
between the density of star formation and the gas den- 
sity and, when including the e/cfit factor, reproduce well 
the observational results, except possibly for the densest 

1 It must be kept in mind that all these SFR values apply to giant 
molecular clouds. Values inferred for entire galaxies, including the 
Milky Way, are substantially lower, as they include diffuse atomic 
or molecular gas, overestimating the amount of gas counted as 
star-forming gas (e.g. Heiderman et al. 2010). 

2 This applies to eqn. ((3)1 as well, whereas the true PN relation 
corresponds to e/cj>t = 1.0 and 0.5 in the hydro and MHD case, re- 
spectively, and should be moved upward accordingly on the figure. 
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Fig. 1. — SFR^j as a function of a„ 



t vlr for various cloud parameters and M = 16 (solid), 9 (dot) and 4 (dash), as predicted by various 

theories. Top panels: cloud size and mass. Second row: Padoan & Nordlund (2011)(left) and the corresponding "multi-freefall PN" 
(right). Third row: Krumholz & McKee (2005) (left) and "multi-freefall KM" (right); Fourth row: Hennebelle & Chabrier complete theory 
(eqn. (0;left) and simplified "multi-freefall" theory (eqn. (|8]|;right) for y cu t = 0.1. 



clumps, where the SFR's are about a factor ~ 5 larger. 
The large spread of the data precludes a clear distinction 
between the theories at this stage, apart from the KM 
one which clearly lies well below most of the data points. 

Interestingly, both the exact HC and "multi free-fall" 
calculations predict a drastic drop in the SFR below 
S c w 110-120 M Q pc~ 2 and a change of slope in the 
oc T,^ relation, with N « 4.8 and iV w 1.6, similar 
to the Kennicutt-Schmidt relation (Nks ~ 1-4), respec- 
tively below and above E c , in very good agreement with 
the observations (Heiderman et al. 2010). This density 
corresponds to a visual extinction Ay ~ 6 (Ak ~ 1). A 
similar density-threshold for significant dense core popu- 
lation has been identified in several surveys (Onishi et al. 
1998, Johnstone et al. 2000, Kirk et al. 2006, Enoch et 
al. 2007, Lada et al. 2010, Andre et al. 2010). Various 
authors (e.g. Johnstone et al. 2004, Kirk et al. 2006, 
Heiderman et al. 2010) have suggested that the origin 



of such a density threshold is related to magnetic fields, 
which cannot support the gas against gravitational col- 
lapse above some density. The present calculations, how- 
ever, show that there is no need to invoke magnetic field 
support or MHD shock conditions to get such a threshold 
although, as mentioned in the appendix, magnetic fields 
may contribute dynamically by reducing the value of 4>t- 
The threshold simply stems from the fact that at the 
corresponding density, the size of the clumps becomes 
comparable to the Jeans length and thus the amount of 
gas appropriate to form stars drops drastically (see HC09 
Fig. 8). The relative similarity between the various 
"multi free- fall" predictions clearly indicates that, be- 
sides the e/<pt factor, the key physical quantity which de- 
termines the SFR is the density- dependence of the freefall 
time of the collapsing overdense regions induced by tur- 
bulence. An alternative possibility, as recently suggested 
by Krumholz et al. (2011) is to assume that the SFR 
is simply a constant factor times the cloud's or galaxy's 
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Fig. 2. — Comparison of the SFR per unit area, as a function of gas surface density, E 9 , as predicted by the various theories 
illustrated in Fig. 1, with the observational determinations of Hcidcrman ct al. (2011) for massive clumps (triangles) and molecular clouds 
(diamonds+squares). The stars show the data of Lada et al. (2010). The four solid lines correspond to four cloud sizes, namely L c = 1, 4, 
10 and 40 pc (left to right). 



volume density over mean free-fall time, although a clear 
physical explanation is lacking at this stage. 

4. CONCLUSION 

We have included the time dependence in our analyt- 
ical theory of the IMF to determine the SFR. The the- 
ory, based on a gravoturbulent picture of star formation, 
yields SFR values in good agreement with various ob- 
servational determinations in Galactic molecular clouds. 
Moreover, it naturally predicts a density threshold to get 
significant star formation and yields a dependence of the 
SFR upon gas surface density in very good agreement 
with the observationally inferred values, with an abrupt 
change of slope around the threshold. Such a threshold 
naturally emerges from our theory, without arbitrarily 
introducing a critical density. 

A crucial point at the heart of the present (both com- 
plete and simplified (eq. (8))) approach is that, in con- 
trast to previous theories, the SFR is not character- 



ized by a single, characteristic dynamical time in the 
cloud, but instead involves a density- dependent collaps- 
ing time for each turbulence-induced gravity-dominated 
overdense region in the cloud. Therefore, in opposite to 
conclusions based on previous SFR theories, the present 
results show that a SFR determined by turbulence- 
induced density fluctuations at the early stages of star 
formation provides quite a consistent picture of star for- 
mation in Milky Way molecular clouds. 
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APPENDIX: THE CROSSING TIME 
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Our estimate for the crossing time is similar to the 
estimate of Krumholz & McKee (2005). The crossing 
time of a structure of scale R is r ct (R) = 2R/V ct . At 
large scales, V ct ~ V imS} while at small scales, below 
the sonic length, V ct ~ C s . The typical time tr within 
which the density field is significantly modified at scale 
R, implying that a new set of fluctuations, statistically 
independent of the former one, has set up is tr = a ct T ct , 
with a ct a dimensionless coefficient of the order of a few. 
In the Hennebelle-Chabrier theory, we select the pieces 
of gas which are self-gravitating. At large scales, this 
implies a g GM/R > V^ ms , where a g is a dimensionless 
coefficient (a g = 3/5 for a uniform density cloud), while 
a similar expression holds below the sonic length. This 



yields 



tr 



2a ct R 

^rms 



where 



2a r 



3tt 



/ 24 
TT 2 a g 



T ff = <t>tTff = 4>tT f} y/p /p, 

(9) 



T ff~\l 32Gp 1S ^ ne free-fall time of a bound region 

of density p, and <pt ~ 3, yielding e/(f> t « 0.1-0.2. 

Note that this estimate assumes that it takes about 
one crossing time to rejuvenate a self-gravitating struc- 
ture. However, it may happen, in particular in magne- 
tized flows, that all perturbations do not collapse eventu- 
ally (e.g. Hennebelle & Perault 2000), further increasing 
4>t by a factor of a few. 
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